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■ We apply the effective field theory definition of Minimal Flavour Violation 

^ ! (MFV) to the MSSM. We exphcitly show how, by this definition, the 

new sources of flavour and CP violation present in the MSSM become 
functions of the SM Yukawa couplings, and cannot be simply set to zero, 
r~| , as is common wisdom in phenomenological MSSM studies that assume 

MFV. 

We apply our approach to the MSSM AB = 2 Hamiltonian at low tan /?. 
The limit of MFV amounts to a striking increase in the predictivity of 
5^ I the model. In particular, SUSY corrections to meson-antimeson mass 

differences AM^^s always found to be positive with respect to the 
SM prediction. This feature is due to an interesting interplay between 
chargino and gluino box diagrams (the dominant contributions) in the 
different mass regimes one can consider. 

Finally, we point out that, due to the presence of gluinos, the MFV MSSM 
does not belong - even at low tan/3 - to the class of models with the 
so-called 'constrained' MFV (CMFV), in which only the SM operator 
{V - A) (g) {V - A) contributes to AM^^^. Consequently, for the MSSM 
and in the general case of MFV, one should not use the Universal Unitarity 
Triangle (UUT), relevant for CMFV models, but a MFV-UT constructed 
from P^Ks ™d \ Vub\ or 7 from tree-level decays. In particular, with the 
measured value of P^Ks^ MFV implies a testable correlation between \ Vub\ 
and 7. With the present high value of \Vub\, MFV favours 7 > 80°. 
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1. Introduction 

The analysis of accelerator data over the last two decades established the success of 
the Standard Model (SM) pattern of flavour and CP violation. If the New Physics 
(NP) introduced to stabilize the electroweak (EW) breaking scale had a sensibly 
different pattern for such violations, it is natural to expect that it would have been 
already visible in the considerable amount of precise data on flavour changing neutral 
current (FCNC) processes available today. The coming years will show whether this 
picture is altered by new data, in particular CP-violation in the B^-system and rare 
K-decays, where large non-CKM sources of flavour and CP-violation are still possible 
within the most popular extensions of the SM. 
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On the theoretical side, the success of the CKM description of FCNC processes 
triggered the idea that the dynamics responsible for the peculiar form of the SM 
Yukawa couplings may be relevant only at energy scales much higher than the one 
typically introduced to stabilize the EW breaking. Such high-energy dynamics would 
then generate only the Yukawa couplings present in the SM and no additional flavour 
violating structures. From the low-energy point of view, the SM Yukawa couplings 
arc then the only 'building blocks' regulating the amount of FCNC and CP violating 
processes, and their form is then raised as a "symmetry requirement" for the flavour 
sector of any candidate extension of the SM at the EW scale [1]. 

The above mentioned idea is known in the literature as Minimal Flavour Violation. 
One could say that this idea in the quark sector has become the more precise, the more 
experimental data tended toward it. A phenomenological definition of MFV, that uses 
(the explicit occurrence of) the CKM matrix as the only source of flavour violation 
and restricts the set of relevant operators in the low-energy effective Hamiltonian to 
the SM ones, has been introduced in [2]. It implies a set of very special relations [2, 3] 
among observables in the flavour sector, that have been extensively tested in the recent 
years. In particular, the unexpected agreement of the so-called Universal Unitarity 
Triangle [2] with the available data [4] has brought MFV to the fore, raising the 
question how to implement it in NP models, whose flavour sector is a priori unrelated 
to the SM one. 

Already at this stage, we would like to emphasize that, while pragmatic and phe- 
nomenologically useful, the definition of MFV introduced in [2], to be called 'con- 
strained MFV (CMFV) [5] in what follows, is not as general as the one in [1], and 
the difference between the two approaches will emerge from our discussion. For this 
reason, in the present paper we will use the general definition of [1] and only at the 
end we will investigate under which assumptions the limit of CMFV can be reached 
in the specific framework of the MSSM. 

Stated loosely, the basic idea to provide a model independent definition of MFV is 
as follows. The only low-energy remnant structures responsible for flavour violation 
are the SM Yukawa couplings, the single ones 'required' at present by experiments. 
Then, the flavour sector of every extension of the SM at the EW scale should be 
minimal flavour violating if its flavour violating 'building blocks' are exclusively the 
SM Yukawa couplings. This idea has been formulated rigorously by the authors of 

[!]• 

On the phenomenological side, the idea of MFV has very often been advocated 
to better constrain models, whose predictivity is spoiled by the large number of 
parameters, as is notably the case for the Minimal Supersymmetric Standard Model 
(MSSM). Focusing on the latter, many studies do already exist in the literature, where 
the MFV paradigm is explicitly advocated. However, many (most, actually) of such 
studies appeared before the "effective theory" deflnition by [1] and use assumptions 
often not complying with such definition. We stress that the latter is the only one 
that can be unambiguously applied in extensions of the SM, since MFV is defined 
through the formal transformation properties of the SM Yukawa couplings and can 
subsequently be applied to any new source of flavour violation. 
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The aim of this paper is twofold. First, we reconsider the quark flavour sector of the 
MSSM and carefully discuss its MFV limit. We explicitly show how the new sources 
of flavour and CP violation present in the MSSM become in this limit functions of 
the SM Yukawa couplings, and cannot be simply dropped in MFV, as often assumed 
in the literature. A simplistic but intuitive picture is that the off-diagonal entries in 
the soft terms - the genuinely new sources of flavour violation in the MSSM - are not 
zero in MFV, but instead 'CKM-like'. 

We then apply our approach to the specific case of the AB = 2 Hamiltonian in the 
MSSM at low tan /?. Such Hamiltonian, responsible for Bg^d meson mixings, provides 

a concrete and phenomcnologically interesting benchmark for the approach itself. 
Explicit implementation of MFV in the MSSM leads to a striking improvement of 
its predictivity. This is obvious if one thinks that the sector introducing the largest 
number of new parameters is notably the soft sector. The latter is now entirely con- 
strained to be proportional to appropriate combinations of the SM Yukawa couplings, 
so that the main unknowns turn out to be the (real) proportionality factors ('MFV 
parameters'), amounting to 12 independent dimensionless parameters. Furthermore, 
one has to fix some real mass scale parameters ('SUSY scales'): the fi parameter, a 
squark mass scale fn and trilinear coupling A, gaugino masses Mi, M2 and Mg and the 
two Higgs soft terms rn//^^ and ttihj- Hence, the mass scales relevant to the AB = 2 
case are in total 8, but we will see that only a subset of them affects non-trivially the 
calculations. Finally, also tan (3 is of course a parameter, but we set it to reference 
(small) values, whose choice does not affect our main findings. 

Concerning the new sources of flavour violation specific to the AB = 2 case, we note 
that contributions from boxes featuring gluinos and neutralinos, usually assumed not 
to enter MFV calculations [6, 7, 8, 9, 10], do actually contribute. The flavour violating 
structures in their couplings are proportional, as mentioned above, to appropriate 
combinations of the SM Yukawa couplings, the combinations being fixed by the very 
definition of MFV [1]. 

The second aim of our paper is a detailed numerical study of the AB = 2 Hamil- 
tonian in the MFV MSSM at low tan/3. In this study, mass scale parameters can 
be fixed to reference values covering all the physically interesting mass scenarios that 
low-scale SUSY could have. Concretely, we fixed to a small (200 GeV) or interme- 
diate (500 GeV) or large value (1000 GeV), covering both signs. For each choice of /j, 
we then chose the squark mass scale to four benchmark values in the range 100 -r- 1000 
GeV and so forth for the other parameters. We considered a total of 48 scenarios. 
Then, the 12 MFV parameters which govern proportionality to the Yukawa matrices 
are left free to float within reasonable intervals. 

Now, for every mass scenario considered, a random scan of the MFV parameters 
allows to generate a range of predictions for the SUSY corrections to the SM meson- 
antimeson mass differences AMg^d- The predicted corrections display a number of 
remarkable features 

i) For each of the mass scenarios considered, corrections turn out to be always 
positive and to float within a relatively narrow range of values when varying 
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MFV parameters. Specifically, in the case of Bg — Bg mixing, corrections are 
in the range 

AM^USY ^ +(0 ^ 2) ps~^ Given the still lar gc error associated 
with the computation of the matrix elements entering the AM^ determination, 
these corrections, at present, are however not large enough to distinguish the 
SM from the MFV MSSM at low tan (3. 

ii) The positiveness of the sum of the SUSY contributions turns out to be caused by 
the interplay between the two dominant of them, namely chargino and gluino 
boxes. This interplay is mainly dictated by the relative importance of the 
H parameter with respect to the other SUSY scales. This can intuitively be 
understood by observing that, if /x is small, it governs the chargino lightest 
mass eigenvalue, whereas large values of /i increase the importance of scalar 
operator contributions coming from gluino boxes. 

iii) By analyzing the single box contributions, we identify four main scenarios for 

the interplay between chargino and gluino contributions. Such scenarios are 
ruled basically by the magnitude of /x and by that of the squark mass scale ra. 
Variation of the other SUSY scales plays only a marginal role in the qualitative 
picture that emerges. 

iv) Since we restrict our analysis to low tan /3, a naive expectation would be that 
most of the contributions be proportional to the SM left-left current operator, 
since the down-quark Yukawa matrix should be negligible. We find departures 
from this picture, arising when ji is not small in magnitude, and due to gluino 
contributions. Responsible for these departures are, in particular, the LR and 
RR submatrices of the down-squark mass matrix. If the down Yukawa is set 
to zero, these submatrices arc respectively zero or proportional to the identity 
matrix. On the other hand, when the down Yukawa is kept, they give rise to 
the bulk of contributions from operators other than the SM one. 

The rest of the paper is organized as follows. In Section 2 we recall the effective 
theory definition of MFV, in the formalism of the MSSM. Then, in Section 3, we apply 
such definition to the flavour sector of the MSSM, in particular to the soft SUSY 
breaking terms, by discussing their MFV relations to the SM Yukawa couplings. In 
Section 4 we then focus on the AS = 2 Hamiltonian in the MSSM at low tan (5. 
We collect here all the basic formulae and discuss the steps needed to evaluate their 
MFV limit in the light of the procedure described in the previous sections. Section 
5 presents our numerical strategy to explore the MFV MSSM predictions for AM^ 
and a detailed accout of our main findings, in particular the features outlined in the 
above points i) to iv). In Section 6, we then elaborate on our findings, by describing 
additional numerical studies performed to clarify the issues emerged, like the role of 
tan (3. Section 7 is devoted to various considerations on the topic of MFV, triggered 
by what we learned from the study carried out in the previous sections. One of such 
reflections concerns the definition of the Universal Unitarity Triangle, which turns 
out not to be a construction always valid in MFV. To this point we devote Section 8. 
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Finally, Section 9 presents our conclusions and outlook. In the Appendix we collect 
the complete list of Wilson coefficients for the = 2 Hamiltonian in the MSSM at 
low tan 13. 

2. Minimal Flavour Violation: effective theory definition 

In a top-down approach, the question how to build up concretely the MFV hypothesis 
in a given NP framework translates into how to define MFV in the presence of new 
flavour violating interactions, a priori unrelated to the SM ones. To understand this 
point, it is first useful to remind the structure of fiavour breaking in the SM. The 
SM flavour symmetry group and its breaking have been first elucidated in [11, 12]. 
Responsible for such breaking are the SM Yukawa couplings, and their transformation 
properties under the flavour group can be identified by requiring (formal) invariance of 
the Yukawa interactions. Flavour violation is recovered as the spurion Yukawa "fields" 
assume their background values. MFV then demands the Yukawa background values 
to be the only structures generating the observed fiavour (and CP) violation. This 
definition, which has the advantage to hold model-independently, has been introduced 
by D'Ambrosio et al. [1]. 

Following this approach, in the context of a given NP model, every new flavour 
violating "coupling" can be classified according to its transformation properties under 
the SM fiavour group, and - if MFV holds - rewritten in terms of combinations of 
the SM Yukawa couplings transforming in the same way. 

Such procedure has been detailed in [1]. Wc now restate it briefly for the MSSM 
with i?-parity, which is our case of interest. The relevant quantity is the superpotential 
W, which reads 

W = eij (vi'HlQ'W' + Yj'^Hi.Q'W'' + Yi' W^L'^ E' + ixW^Hi) , (1) 

with the matter superfields Q, U ^ D, L, E (containing SM fermions) and (con- 
taining the Higgs doublets). Here /, J and i,j denote flavour and SU{2)l indices, 
respectively. The notation and conventions comply entirely with [13]. 

Out of the largest possible group Gi? of fleld redefinitions that commutes with the 
gauge group [11], 

Gf = [SU{3) ® U{l)f = (g) [SU{3) ® U{l)]p , (2) 

F=Q,U,D,L,E 

the Yukawa interactions in the superpotential (1) break the fiavour group [SU{3)]^ (g) 
U{1)e [14, 1]. The flavour symmetry can formally be recovered in eq. (1) by treating 
Y^^d^f, as spurious and requiring them to have indices transforming under [SU{3)]^ as 
follows 

with the subscript Q,U, D, L, E referring to an index that transforms as the corre- 
sponding representation under SU{3)q^u,d,l,e, respectively, and as a singlet under all 
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the other group factors. Note that the superfields U,D,E are left-handed but must 

describe right-handed particles. As a consequence their component fields are defined 
with a charge conjugation operation and they transform as 3 representations under 
SU{3)u,D,E, respectively. 

Using the [^[/(S)]^ symmetry, the fermion superfields can be suitably shifted to 
have Yukawa couplings in the form 

Y^ = K^Yu, Yd = Yd, Ye = t, (4) 

with the Y diagonal matrices and K the CKM matrix. This form is not the usual one, 
since quark mass matrices are not simultaneously diagonal. However, it is very useful 
when ranking different flavour changing effects, since the top Yukawa (the dominant 
one) displays explicit proportionality to the CKM matrix. 

For low tan/3 = [Hy) / [Hd) , all FCNC effects are dominantly described by one 
single off-diagonal structure [1] 

(YuYi)ij^\lK^iK*^^, i^j, 
I , i = j, 

with At = {Yu)33. Note in fact that higher powers of Y^Yu can be rewritten in terms 
of YuYu times an appropriate power of Xf. Subleading effects on the r.h.s. of eq. (5) 
arc suppressed by powers of nic/mf. 

The main observation [1] is now that, if MFV holds, soft SUSY-breaking terms are 
related to the SM Yukawa couplings (4) and the explicit relations can be constructed 
by just using the formal transformation properties of (5) under the flavour group. 
Such derivation will be presented in section 3. In section 4 we will subsequently use 
the calculated SUSY parameters to evaluate the MSSM contributions to the AB = 2 
effective Hamiltonian, to which we will restrict the rest of the analysis. The latter 
will allow us to quantitatively assess the features of the SUSY contributions in the 
MFV limit in a benchmark case like meson oscillations. 

In this respect, wc anticipate that even if (AFc)i/ in cq. (5) provides the dominant 
FC mechanism, a detailed study shows that effects proportional to the down-quark 
Yukawa cannot actually be neglected. The latter does not provide by itself an addi- 
tional FC mechanism - as one can see from eq. (4) - but still its effects can correct the 
magnitude of those provided by (AFc)ij- This can be understood by simply looking 
at the LR entries of the down-squark mass matrix (see eq. (9) below). The latter 
are proportional to fxYd, and if f^ is not small, the corresponding terms cannot be 
dropped even if <C 1^. 

Focusing again on the AB = 2 Hamiltonian, this mechanism also regulates the rel- 
ative importance of contributions to operators beyond the SM left-left vector operator 
Qi. This can also be naively understood from the down-squark mass matrix, once 
its MFV limit is performed. In this Umit, in fact, the LR and RR sectors are zero 
or respectively proportional to the identity matrix, if Y^ ^0. As a consequence one 
would expect the SUSY contributions to the AB = 2 Hamiltonian be in the Wilson 
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coefHcient of Qi only. This picture is largely modified or completely spoiled when one 

includes Y^, depending on the mass scenario chosen for the SUSY parameters. The 
main actor in this respect is again the /j, parameter, for the reason outlined above. 



3. MFV relations of MSSM parameters to SM Yukawa 
couplings 

We now discuss how the above picture concretely applies to the MSSM with i?-parity 
and softly broken SUSY. In this model, the part of the Lagrangian responsible for 
flavour violation reads 



ni.v. 
■'-MSSM 



with W the superpotential of eq. (1) and JO-l'^fi given as 



(6) 



/•f.v. 

"■'-soft — 



- 'e,^{Ai'Wjj^{uir+Ai'h^Ji{dir 



(7) 



i.e. the usual soft Lagrangian with omitted gaugino mass terms. 

The MSSM Lagrangian, with the above superpotential and soft pieces, gives rise 
to the following 6x6 squark mass matrices 



Ml 



\Y^Yl + (mlf 



4M, 



w) 



-ii*^Y 



^tan/3A 



-H^Y^-^is^npAi 



^Y^Y^ + ml + \ cos 2/3M|4l 



A^^tan/jyJ + ^At 



;x*^tan/3y, + ^A, 



ml - ^M|4l 



If MFV holds, the new sources of flavour violation present in the soft terms must 
be related to SM Yukawa couplings. To this end, they can again be formally treated 
as spurion fields, with indices transforming under the flavour group as follows: 



Recalling Yukawa transformations (3), one can then write the following MFV relations 

ml = m'^ (^a2l + b^Y^Yu^ , 
ml = (agl + fte^Jyd) , 

= A(^a4Y^ + b7YdY}Yu) , 

Ad = A(^a5Yd + bsY,Y^Yd) , (11) 

where, in the first line, we have reported the transpose of m^, since it is the latter to 
appear in the squark mass matrix in the usual conventions [13]. The Oj, 6j coefficients 
are real proportionality factors, whose allowed range of values will be studied in 
section 5 below. The overall mass scales fn and A fix the order of magnitude of the 
respective soft terms, when the coefficient multiplying them is of 0(1). Expansions 
(11) are accurate to all orders in At and = (1^)33. Subleading effects are suppressed 
by powers of mc/mt and/or mg/mi,. 

Starting from the mass matrices (8)-(9), and expressing the soft terms according to 
expansions (11), it is then customary to perform a superficld redefinition leading to 
diagonal mass matrices for the quarks. The unitary matrices adopted are the same 
as in the SM: here one defines shifts 

UL Vq^ul , dL VQidL , 

UR VuUR , dn Vodn , (12) 

and the fields on the r.h.s. form the CKM basis. In the MSSM, the same shifts are 
carried out at the superfield level and lead to the so-called super-CKM basis. After 
performing such transformations, one gets diagonal Yukawa matrices 1^, 1^, and can 
use relations 

mu = -^Yu, md = ~--^Yd, (13) 

in eqs. (8)-(9) to display explicit dependence on the quark mass matrices ifiu^d- 

In the super-CKM basis, the matrices mg ^ ^ and A^^d have still off-diagonal en- 
tries.^ Then, in order to have the (hermitian) mass matrices in eqs. (8)- (9) in diagonal 
form, one needs a second redefinition, performed on the up- and down-squark fields, 
respectively. Such redefinition leads from the basis u, d, to the mass eigenstate basis 
U, D, which in the conventions of [13] reads 

Ui = iZu)ijUj , di = {Z*o)ijDj , (14) 



^Note that 63 and 64 must be equal, due to the hermiticity of mq. 

^Such entries are responsible in general for genuinely supcrsymmetric flavour violation in the mass 
matrices (8)- (9). In our case, as we said, soft term are instead fixed by the MFV expansions (11). 
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with the index i = 1,2,3 for ul, cLl and z = 4, 5, 6 for ur, dji. 

Using transformation (14), the down-squark mass term in the Lagrangian is diag- 
onalized according to 

FMjd* = D'^[z],MjZD]D* = D'^MlD* , = diag{Ml^, M^J , (15) 

with M| given in eq. (9) and d {D) a column vector built out of the dj (Di). An 
entirely analogous equation holds for the case of up-squarks. 

In practical calculations, flavour violation is driven either by non-diagonal squark 
propagator matrices, when working in the d, u basis, or by a Zd,u niatrix appearing in 
vertices with a squark leg (gluino-quark-down squark, neutralino-quark-down squark 
and chargino-quark-up squark, in our case). In the former case, one usually adopts 
an expansion in off-diagonal "mass insertions" and stops to the first non-trivial order, 
in the well-know Mass Insertion Approximation (MIA) [15, 16]. The MIA provides a 
very useful tool to make flavour violation mediated by soft SUSY breaJiing terms most 
transparent and manageable, since it 'linearizes' the mechanism of flavour violation, 
but it is an approximation. 

In our MFV formulae, we will instead stick to the mass eigenstate basis, i.e. to 
the exact calculation. Among the M| entries in eq. (9), the soft terms will be 
related through proportionality factors to the SM Yukawa couplings, according to 
the expansions (11). Then the Mj mass matrix turns out to depend only on fi (which 
must be real), on tan/J, on the two squark scale factors A,rn (see eq. (11)) and 
on the proportionality factors. Upon rotations of the squark states from the super- 
CKM basis to the mass cigcnbasis, eqs. (14), the pattern of flavour violation is then 
transferred from the non-diagonality of the mass matrices, to the off-diagonal entries 
of the matrices Zj^^u, entering quark-squark interactions with gluinos and neutralinos. 

Let us show with a simple example how flavour violation in the Zj)fj becomes 
'CKM-like', after MFV expansions are imposed. Let us consider the down-squark 
mass matrix of eq. (9), with soft SUSY parameters given according to the MFV 
expansions in eq. (11). Adopting the approximation Y^i 0, one can drop all 
the corresponding terms in such expansions. One can then perform the super-CKM 
rotation on the squark flelds to have the up Yukawa diagonal. The down-squark mass 
matrix assumes, in this basis and under these assumptions, the following form 

^ / m\ar\ + hriK'^n^f) - ^{Ml + 2M2.)1 \ 



whence the unitary transformation (14) leading to the mass eigenbasis for the down- 
squarks is obviously 
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As one can see, ofF-diagonal entries in are not zero, but CKM-like (in this sim- 
ple case only in the LL sector). They will appear in the couplings of gluinos and 
neutralinos with quarks and squarks. 

When one includes in the diagonalization the effects of the down- Yukawa matrix, 
the diagonalization becomes more involved. However, flavour violation in the ^d,!7 
is still encoded in their dependence on the SM Yukawa couplings. 

This observation allows us to comment on a conventional assumption present in 
most of the calculations performed in the MFV MSSM to date. This assumption 
amounts to dropping altogether flavour violating entries in the Zjj^u, due to the com- 
mon wisdom that, if MFV holds, flavour violation can come only from couplings 
explicitly displaying proportionality to the CKM matrix. The effective theory defini- 
tion of MFV [1] implies that the correct approach is instead to think the off-diagonal 
terms in the Zd,u as being not zero, but instead dictated by the SM Yukawa cou- 
plings. In the simple example above, this dependence reconstructs in (the LL sector 
of) Zf) directly the CKM matrix. 

The bottom line is that, in the MFV MSSM, one has to diagonalize squark matri- 
ces after imposing expansions (11), so that the diagonalization matrices Zd,u bear 
dependence on such expansions and then use the Zu^jj in all vertices where they are 
required. This point has already been stated in [1]. 

For the sake of completeness, we also report here the chargino, ncutralino and 
charged Higgs mass matrices, since these particles will enter our subsequent calcula- 
tions. The notation used is again that of [13]. Charginos axe two Dirac fermions xi,2 
whose mass matrix reads 



M. 



M, 



X2 




(18) 



with Z± unitary matrices, chosen from the requirement < M^^ < M^^. Similarly, 
neutralinos are four Majorana fermions Xi 4; with mass matrix given by 









( Ml 


—evi 




Ms 

evi 


—evi 









\ 




\/2cw 

ev\ 


\p2cw 
—ev2 








— 


\f2syif 




\/2sw 




V 


^x2 


/ 




ev2 


—ev2 




\ V2cw 


\Pisw 



z 



N 



(19) 



with Zn a unitary matrix, whose form is again specified after requiring positiveness 
and ordering for the eigenvalues. 

Finally, one has two physical charged Higgs scalars Hf, with mass 



(20) 



where fn^j^ and 'rn?H^ are soft terms for the corresponding Higgs doublets, given in 
eq. (7). Away from the unitary gauge, one must also include in the calculations the 
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fields, which provide the longitudinal degrees of freedom for the W bosons in the 

unitary gauge. 

When assuming MFV, the gaugino masses Mi_2 are real.^ In fact, if one allows 
non-trivial phases in Mi ^2, they are communicated to the diagonalization matrices 
Zn and Z±, which in turn enter Feynman rules for charginos and neutralinos. One 
would then have new sources of CP violation, not allowed by the MFV hypothesis. 
The same argument applies to the Higgs sector parameter /i. 

4. AB = 2 in the MFV MSSM at low tan/3 

The MFV limit of the MSSM, as described in the previous sections, can now be 
applied to a concrete example, that of the AB = 2 Hamiltonian, which is responsible 
for meson-antimeson oscillations. The latter has recently received renewed theoretical 
interest, in view of the very precise measurement of Bg — Bg oscillations by the CDF 
collaboration [17]. 

The basic ingredient to describe meson-antimeson oscillations is the quantity Mi2^ 
= {M\n^/='^\M), with M = K, Ba,s- Within the MSSM, H^/=^ has the form 

5 3 

= Y,CiQi + ^QQi + H.c. , (21) 

i=l i=l 
with the Qi given, in the case of Bg — Bs mixing, by 

Q2 = {s'Pl b') {s^Pl V) , 

Q3 = (s'Pl V) (s'Pl b') , 

Q4 = i^'PLb'){-s^PRV) , 

Q5 = (^'PLV){-s^PRh^) . (22) 

The operators Qi.2,3 are obtained from Qi,2,3 by the replacement L R. The left- 
and right-handed projectors are defined as Pr^l = (1 =t75)/2 and = '^'^Pr,l'i hj 
are colour indices. In the case of 5^, one should replace s — d in eq. (22). 

Each of the Wilson coefficients in eq. (21) features, for low tan/3, the following 
contributions 

Q = Cf^ + cf^H+ + ^x+x+ + + ° + Cf^' , (23) 

where, for simplicity, we have omitted to specify the flavour indices of the external 
quarks, as in eq. (21). In eq. (23), the first term on the r.h.s. represents contributions 

from the SM boxes. The additional contributions, that need to be considered within 
the MSSM, come respectively from boxes with charged Higgs-up quarks, chargino- 
up squarks, gluino-down squarks, mixed gluino- and neutralino-down squarks, and 

^Mg can be chosen as real without loss of generality [13]. 
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Figure 1: Feynman diagrams describing meson- antimeson oscillations in the MSSM. The 
crossed diagrams (second row) are needed only if the fermion in the loop is a 
Major ana particle. The notation for the various lines is the same as in [13]. 

neutralino-down squarks. The possible Feynman diagrams involved in each case are 
represented in Fig. 1. For the SM, charged Higgs and chargino cases, one has Dirac 
fermions propagating in the diagrams, so that only the boxes in the first row of Fig. 1 
must be considered. The other contributions involve Majorana fermions in the loop, 
so that also crossed boxes (second row of Fig. 1) need to be calculated. 

The complete list of Wilson coefficients for the AB = 2 effective Hamiltonian in 
the MSSM is reported explicitly in the Appendix. Such coefficients are calculated in 
terms of loop functions, depending on the particle masses involved in the loops, and 
couplings, possibly featuring the rotation matrices Zu^d (squarks), Z± (charginos), 
Zn (neutralinos) introduced in the previous Section to define the respective mass 
eigenstates. To evaluate the Wilson coefficients in MFV, the procedure to follow is 
now clear: 

1. Expand the soft terms as in eq. (11) if they transform non-trivially under the 
flavour group, or take them as real if they are singlets; 

2. Plug them into the mass matrices and diagonalize the latter to obtain the mass 
eigenvalues and the rotation matrices defining the eigenbases; 

3. Use the obtained eigenvalues and rotation matrices in the general MSSM for- 
mulae for the Wilson coefficients. 

Now that we have all the ingredients of the calculation, we conclude by listing the 
number of parameters involved in the MFV limit of the AB = 2 MSSM Hamiltonian 
for low tan /?. The expansions of the soft terms in the squark mass matrices, eq. (11), 
involve 12 real proportionality factors, and 2 overall mass scales: a 'generic' squark 
mass m and a 'generic' trilinear mass term A. In addition one has to fix three real 
gaugino mass terms Mi, M2 and Mg and the real fi parameter. Finally, the soft 
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Higgs sector adds 2 more mass scales, namely niH^ and iJiHa- Taking into account 

the requirement of correct EW symmetry breaking, amounting to one constraint, one 
has a total of 12 + 7 parameters.^ Among the latter, actually the dependence of the 
computed Hamiltonian on mn^ ^ t)e trivially 'factored out' and not considered 
in MonteCarlo approaches. In fact, mj/^ ^ enter only Higgs boxes, and the latter 

in turn depend exclusively on ^ and on n, through eq. (20), which fixes the 
physical charged Higgs mass M^±. There is instead no dependence on any of the 

other SUSY scales and on the MFV parameters. In addition, considering that ^ 
are demanded by the soft Lagrangian (7) to be real, but not necessarily positive, it is 
clear that, for every choice of /i, it is always possible to tune ^ (compatibly with 
the EW symmetry breaking constraint) so that M^± assumes any desired value. As 
a consequence, one can trade the parametric dependence on ^ for that on M^± , 
which is then the only SUSY parameter in the Higgs contributions. 

In the next section, we will discuss the MonteCarlo procedure adopted to explore 
the above parameter space. We will see that the relevant quantities to be scanned 
turn out to be only the 12 MFV proportionality factors, so that the predictivity and 
testability of the model end up to be dramatically improved. 

5. MFV MSSM predictions for meson mixings 

We are now ready to study the MFV MSSM AB = 2 amplitude, and its predictions 
for AMs^d, by varying the SUSY scales as well as the MFV proportionality parame- 
ters. To this end, a MonteCarlo approach which generates every (or a subset) of the 
parameters according to fiat distributions within given ranges provides the most sys- 
tematic and unprejudiced tool. Here below we describe in more details our adopted 
strategy. 

5.1. Strategy 

Our reference numerical study was carried out by fixing the mass parameters to 'sce- 
narios', and then, for each scenario, scanning with fiat distributions the 12 parameters 
(aj,6j) ruhng the MFV expansions (11). In addition, we set tan/? = 3.^ 

The choice of the mass parameters was designed to cover, in an exhaustive way, 
all the mass scenarios reasonably allowed for low-energy SUSY. To this end, we have 
started from considering the information on the ranges permitted to SUSY masses, 
that is provided by experiments [18]. On this point we make the following remarks 

• Concerning squark masses and Mg, the most updated bounds ([19, 20] and 
updates thereof) are given in the plane Mq — Mg, where Mq denotes a generic 

*We note that soft terms, expanded according to eq. (11), do actually depend on the product 
between a mass scale and a MFV coefficient, so that the real parametric dependence is on the 
product between the two. Considering this, the above counting is somehow an overcounting. 

^The impact of variations of tan/3 in the range [3, 10] was addressed in a specific set of runs to be 
described below. 
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squark mass (see [19, 20]). The profile of the bound is such that when the gluino 

mass can be small, then the generic squark mass is constrained to be large, and 
viceversa. We have chosen four points in the plane Mq — Mg, namely ^ 

{Mq, Mg) = {(100, 700), (200, 500), (300, 300), (1000, 195)} GeV , (24) 

and used them to fix respectively fn and Mg. We note that m is strictly a repre- 
sentative quantity for the squark mass only when the a-parameter multiplying 
it (sec cq. (11)) is of 0(1). However, the detailed choice of m turns out to play 
a marginal role in our main findings, and the above argument serves only to 
give a reasonable criterion on choosing the pair of parameters fn and Mg. We 
further note that the case in which both fn and Mg are small is 'covered' when 
the a-parameter is small. 

• The generic trilinear coupling A was fixed to the value A = 2fn. ^ This choice, 

when TTT is a representative quantity for squark masses, helps having a not too 
low mass for the lightest Higgs [21]. Again, the full spectrum of deviations from 
this relation is actually covered when scanning the a-parameters multiplying fn 
and A. 

• Constraints on ^ are generically model-dependent. We then considered small, 
intermediate and high values for its magnitude by setting the following possi- 
bilities 

/X = {±200, ±500, ±1000} GeV . (25) 

• Ml and M2, as well as fi, enter chargino and neutralino mass matrices. The 
choice of Mi and M2 in connection with that of fi determines the amount of 
gaugino- and higgsino-like components in their field content. In order to have 
representative cases with respect to the experimental information [18] , we made 
the following choices 

1^1 = 200 ^ (Ml, M2) = {(500, 500), (1000, 1000)} GeV , 

= 500 ^ (Mi,M2) = {(100, 200), (500, 500)} GeV , (26) 
1^1 = 1000 => (Ml, M2) = {(100, 200), (100, 500)} GeV , 

i. e. two possible choices for every of the six /x values listed in eq. (25) . Choices 
(26) translate into values for the masses of the lightest chargino and neutralino, 
which in turn tune the importance of the respective box contributions. We 

^Tho bounds in the {Mq,Mg) plane provided by Refs. [19, 20] are in fact somehow tighter than 
the values chosen in eq. (24). We note however that these experimental bounds are obtained 
assuming a specific mSUGRA scenario. Moreover, in the present study we take the approach of 
preferring smaller masses, in order to address the possibility of large signals in meson mixings. 
As it will emerge from the discussion, even in this approach NP signals in the MFV MSSM are 
however typically found to be within present errors associated with mixings themselves. 

■^For m = 1000 GeV we chose also A = 1000 GeV. 
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note that, in our case, neutralinos have almost no impact on the sum of the 
contributions, even when they are very hght. This apphes also to the mixed 
gluino-neutralino boxes. In this respect, we observe that, from eqs. (24)-(26), 
there are unphysical cases among our considered scenarios in which the light- 
est SUSY particle (LSP) is not a neutralino. However, given the mentioned 
marginal impact of neutralino masses on our main findings, one can always 
lower the value of Mi in order to have a neutralino as the LSP. The values 
chosen in eq. (26) are intended to ascertain that the impact of the choice of the 
neutralino masses be in fact minimal. 

Concerning charginos, their contribution is regulated by the lightest between 
M2 and |//| (see eq. (18)), the detailed choice of the other parameters playing 
basically no role on the main findings we will discuss. 

• The remaining two parameters ^ enter exclusively Higgs boxes. As we also 
remarked at the end of Section 4, the calculation of the latter can be 'factored 
out', since they depend only on mn^^ A*' through relation (20), and 

on no other SUSY scale. With reference still to the discussion at the end of 
Section 4, it is also clear that the single relevant SUSY scale introduced by the 
Higgs sector is the physical charged Higgs mass M^±, eq. (20), and not ruH^ ^ 
separately. 

The dependence of Higgs contributions on variations of and on tan /3 will 
be studied in a separate Section below. Here we mention that such contributions 
are positive for every allowed value of M^± if tan /3 < 7, and even for tan /3 = 10, 
they reach (small) negative values only with very light M^± . As a consequence, 
their impact for low tan/3 is just an overall (positive) shift of the sum of the 
other contributions. 

This completes the discussion on the choice of the SUSY scales in our main analysis. 
The mass scenarios explored amount to 48. Taking into account the various remarks 
made above on every specific subset of the parameters, we believe that such analysis 
covers extensively all the interesting combinations in the SUSY parameter space. 

For each of the above scenarios, we then scanned the MFV parameters 
suming (uncorrelated) flat distributions according to (see also footnote 1) 

0.25 < ai,2,3 < 1 , - 1 < Ks, bi,...,8} < 1 . (27) 

The lower bound in ai^2,3 was not chosen to be zero, in order not to have to discard 
most of the resulting squark matrices because of a negative lowest eigenvalue. 

We finally note that, in our analysis, we do not include other FCNC constraints 
which could in principle play a role for low tan /3, in particular h 57. In this respect 
we observe that, as already mentioned in the introduction, NP corrections to meson 
mixings within the MFV MSSM are typically within present errors and the inclusion 
of additional constraints can only further suppress NP signals. In addition, as again 
mentioned in the introduction, our analysis will lead to the identification of mass 
regimes, ruled by the interplay between chargino and gluino contributions, with the 
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Higgs contributions discussed separately. Since in b ^ the main role is played by 
chargino and Higgs contributions (see e.g. [22]), it is clear that the 6 — ^ 57 constraint 
would not exclude any of the above regimes. 

5.2. Results 

We now discuss our results. The latter were all obtained using the MonteCarlo 
strategy outlined in the previous discussion. We have however also verified the specific 
findings with alternative runs, designed to uncover possible loopholes. We will refer 
to them in due course. 

Our phenomenological analysis starts from the calculated meson-antimeson oscil- 
lation amplitude 

^ {M\Ht/='m , (28) 

with M = K, Ba^s- Wilson coefficients, evaluated at the matching scale, are subse- 
quently run to the rrih pole mass or to 2 GeV, the scales at which the effective matrix 
elements arc evaluated on the lattice [23, 24] in the B^ g and K case, respectively (see 
also [25, 26, 27, 28, 29]). In the running, we used NLO formulae from [30, 31], with 
the matching scale chosen to be at 350 GeV, as a compromise between all the SUSY 
scales entering the calculation.^ 

Then, from twice the amplitude (28), one can calculate the experimentally mea- 
sured mass differences by taking the absolute value or respectively the real part in 
the Bd^s or K cases [32]. In the present study, we restrict to AMg^d- 

As a first check, we have verified that the phase of meson-antimeson oscillations 
in the MFV MSSM be aligned with the SM one. This feature is shown in Fig. 2 for 
the AMg case. In the Figure, we have chosen a specific mass scenario for the SUSY 
scales and scanned the MFV parameters Oj, 6,, obtaining a distribution of values for 
Aii2- As expected, values are aligned along a line with slope tan(arg A^i2). 

In MFV, the phase of meson-antimeson oscillations is by definition not a good 
observable to search for NP effects. However, the same does not apply to the mass 
differences. As a matter of fact, by studying the latter, we found a number of inter- 
esting and sometimes surprising features, which we now discuss. 

i. NP contributions are positive 

A first surprising fact emerges by studying the sum of the SUSY contributions to the 
meson mass differences. As already mentioned above, for every of the mass scenarios 
considered, we have randomly generated the MFV parameters Oj, bi. The obtained 
distribution of values in M12 translates into a corresponding distribution for the 
meson mass differences. As an example, one can look again at Fig. 2, displaying M12 
for the i?s-meson. In this case, Im(A^i2) — and an excellent estimate of AM^^ is 
provided by the projection of the distribution along the Re(A^i2). In the left panels 

"^Variations around this value have basically no effect on the results. 
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99.10 % positive 




Figure 2: Lego plot showing the ahgnment of the phase of SUSY contributions to the SM 
phase in the MFV MSSM M12 for the case of Bg- In this example SUSY scales 
are (GeV): n = 1000, m = 200, Mg = 500, Mi = 100, M2 = 500. In the legend, 
AMs is calculated from the absolute value formula, keeping the sign of the real 
part. The percentage gives the integrated number of hits for which AAlf^ > 0. 

of Figs. 3 to 6 we show how the distribution in AM^^ looks like in four representative 
scenarios.^ 

As one can immediately realize, almost the totality of points features AM^^ > 0, 
i.e. SUSY contributions to AMs in the MFV MSSM at low tan /3 are positive. We 
explicitly mention that in Fig. 2, as well as in Figures 3-6 below, charged Higgs 
boxes are not included. Their contribution amounts to a further positive shift of the 
distribution representing the sum of the other contributions, since Higgs boxes do not 
depend on the MFV parameters , 6j . 

In the following discussion, we will get further insight on the positiveness of the 
SUSY corrections to meson mixings, by analyzing the separate SUSY contributions 
which sum up to give AM^^ . Their interplay and a number of additional checks turn 
out to provide as many arguments in support of the above statement. 



"in the plots, the number of 'events' obtained after scanning at, bi is set to 1000. We have verified 
that the distributions are left quahtatively unchanged when considering subsets of these 1000 
points and when changing the binning, so that 1000 is a statistically significant number. 
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Figure 3: Distribution of values for AMf ^ in the MFV MSSM: sum of the contributions (left 
panel) and separate SUSY contributions (right panel). The distribution results 
from scanning the MFV parameters Qi, bi, after choosing SUSY scales as (GeV): 
^j. = 200, m = 300, Mg = 300, Mi = 500, Ah = 500. In the plot and in the legend, 
AMg is calculated from the absolute value formula, keeping the sign of the real 
part. The percentage gives the integrated number of hits for which AMf ^ > 0. 
[See also text, mass regime (A).] 



ii. Mass regimes 

The positiveness of the sum of SUSY contributions holds true, irrespective of the mass 
scenario chosen. As a matter of fact, the 48 scenarios we have considered turn out to 
be classifiable into 4 main 'mass regimes', each characterized by a definite interplay 
among the various SUSY contributions. The deciding factors are the choice of the 
generic squark mass m and of the magnitude of ^, the rest of the mass parameters, 
as well as the sign of fx, playing only a minor role once these are fixed. Specifically, 
one can choose a 'not large' value for the generic squark scale m (i.e. m < 1000 
GeV in our eq. (24)) and (A) small or (B) |;u| intermediate or (C) \fj,\ large. 
Alternatively, there is one more mass regime (D) when one chooses a large value 
for m, irrespective of the choice of all the other SUSY scales, including fi. The four 
scenarios displayed in Figs. 3 to 6 are representative of such mass regimes, in the 
order (A) to (D). 

iii. Interplay between chargino and gluino boxes 

Let us now have a closer look into the various mass regimes, by discussing, within each 
of them, the main features of the separate contributions. The latter are displayed 
in the right panels of Figs. 3 to 6. As one can see from the figures - and as it 
will emerge from the subsequent discussion - the main actors in determining the 



-19- 




Figure 4: Same as Fig. 3 but for SUSY scales chosen as (GeV): fi = 500, m = 300, Mg = 300, 
Ml = 100, M2 = 200. [See also text, mass regime (B).] 



sum of SUSY contributions to meson oscillations are chargino and gluino boxes, the 
remaining ones playing a minor role. We remind the reader that Higgs contributions 
are not considered in this discussion and not included in Figs. 3 to 6. For low tan /?, 
Higgs contributions trivially amount to a positive shift of the total result. We will 
address this point in a specific section. 

Mass regime (A). Here, the smallness of ji implies a low value for the light- 
est chargino mass eigenstate (see eq. (18)). As a consequence, contributions from 
charginos tend to be large, of the order of 2 ps~^. In addition they are positive, being 
dominated by the SM operator Qi [6]. Gluino contributions are negligible in this 
scenario, and the reason, still related with the smallness of will be clear from the 
discussion in mass regime (C) below. 

A decrease in m only reinforces the chargino dominance. One finds in this case that 
gluino boxes remain negligible, while, for chargino ones, the up-squark mass scale 
gets lower and their contribution is correspondingly increased. As a matter of fact, 
Fig. 3 shows somehow the 'worst' case among the low |^| ones within our studied 
scenarios. In the other cases, the chargino dominance is even more evident and the 
number of points with AM^^ > even closer to 100 %. 

Mass regime (B), with |/i| moderate, is a case of transition, intended to show the 
rate of decrease in importance of chargino contributions with increasing Chargino 
contributions are still dominant as compared to gluino ones. Hence, in this respect, 
the situation is not qualitatively different from regime (A). On the whole, an increase 
in /.i from 200 to 500 GeV corresponds to a decrease in the total signal from ~ 2 to 
~ 1 ps~^, and similarly to regime (A), the total signal drops abruptly beneath the 
peak value. 
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Figure 5: Same as Fig. 3 but for SUSY scales chosen as (GeV): — 1000, m = 300, 
Mg = 300, Ml = 100, M2 = 500. [See also text, mass regime (C).] 



Mass regime (C) occurs then for large In this case, the flavour diagonal 
LR entries in the down-squark mass matrix are large (see eq. (9)) and in connection 
with the flavour mixing induced by LL entries, enhance contributions from scalar 
operators in gluino boxes. Even if the up-squark mass matrix has in principle the same 
structure, the same mechanism turns out to be not efficient in enhancing chargino 
scalar contributions.^'^ The main parametric dependence ruling the above mechanism 
for gluino contributions is then the product fi x tan/3 (see eq. (9)). One should 
also note that the increase in the value of with respect to the previous regimes 
largely suppresses chargino contributions, helping in turn chargino and gluino boxes 
to become of comparable size. 

The right panel of Fig. 5 displays a typical case for mass regime (C), with gluino 
contributions amounting to roughly 30 50 % positive corrections to the chargino 
signal. We note here the occurrence of another interesting mechanism: in cases where 
the squark scale fn and/or the mass Mi are small (as in the example of Fig. 5), gluino- 
neutralino boxes give a negative and relatively important contribution. However, the 
latter is typically outpaced by the positive contribution from pure gluino boxes (plus of 
course that from charginos), with a total signal around 0.5 ps~^. A similar mechanism 

^"One can provide an intuitive argument for this fact as follows: gluino contributions to the various 
operators have the structure ZoZoiyoop function) Zj^Zj^, with Zd defined in eq. (14). Such 
structure holds for every operator. On the other hand, chargino contributions to scalar operators 
have the structure Vi Vi{(loop function) V^Vj^, with Vl(r) the left (right) chargino-up squark-down 
quark vertex coupling (see Appendix), while a similar structure - but with four V^^ vertices - 
holds for the contributions to Q\. Now, since Vl ~ i^u and Vr ~ Yd (see Appendix A. 2), in the 
case of charginos, contributions other than Q\ are always made small by a suppression factor of 
{Yd/Y^f. 
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Figure 6: Same as Fig. 3 but for SUSY scales chosen as (GeV): — 500, m = 1000, 
Mg = 195, Ml = 100, M2 = 200. [See also text, mass regime (D).] 



can already be recognized in regime (B) (see Fig. 4), but in that case is less evident. 
We mention that, within the set of scenarios corresponding to regime (C), we found 
'extreme' cases where scalar contributions from gluino boxes completely overwhelm 
any other contribution. These occur when choosing a very light squark scale, m = 100 
GeV. In these cases, chargino boxes amount to a small positive signal, while gluino- 
neutralino boxes give a contribution which is negative and relatively large. The latter 
is again significantly counterbalanced by the positive, large signal from gluinos and 
the sum of contributions results in a positive, quite spread signal for AM^^ . However, 
such extreme cases correspond to very light squark masses, the lightest down- and 
up-squarks being around 30 and 60 GeV, respectively, which is very unrealistic.^^ 

Mass regime (D) is characterized by a large value for the squark scale m, in our 
case 1000 GeV, and basically unaffected by the choice of any other parameter. In 
this case, the largeness of the squark scale sets to zero gluino-neutralino contributions, 
whose negative skewness had some effect in the previous cases. On the other hand, 
both chargino and gluino distributions (positive) are characterized by a long tail, as 
shown on the right panel of Fig. 6. The respective contributions are comparable in 
size, with the magnitude of those from gluinos growing with The average SUSY 
signal is generically small, ^ 1 ps^^. 



As a further remark, we note that choosing \^\ large, with m small, causes the squark mass 
matrices' determinants to be negative for most of the parameter space in the MFV parameters 
ai,bi, implying in turn an odd number of negative mass squared eigenvalues. The corresponding 
point in the parameter space is then discarded as unphysical. As a matter of fact, for = 1000 
GeV and m = 100 GeV, the ratio of discarded to valid points is ~ 125, but this number drops to 
<; 2 already for m = 200 GeV. 



-22- 



As a last overall remark, we explicitly note that, in all regimes considered, the 
magnitude of contributions coming from the flavour off-diagonal elements in the down- 
squark matrix, entering gluino and neutralino boxes, typically docs not exceed 0.5 
ps^^. Therefore the latter set of contributions is strictly important only when also 
chargino contributions are small, i.e. in regimes (C) and in particular (D). We may 
add that such regimes are phenomenologically relevant after the simple observation 
that the experimental measurement of Bg oscillations [17] and its agreement with the 
SM central value undoubtedly favour small NP corrections with respect to large ones. 
The final word will be provided by a substantial decrease of the lattice error, to the 
level of a few percent. 

iv. Role of charged Higgs boxes 

In the above study, we have completely omitted the inclusion of charged Higgs boxes. 
Recalling the discussion at the end of Section 4, we can now investigate their contri- 
bution separately, as a function of the single new SUSY scale they introduce, namely 
the physical charged Higgs mass M^± . 

In Fig. 7 we report the charged Higgs contribution to AMg as a function of M^± , 
for values of tan/3 between 3 and 10. The hatched area on the left of Mjj± = 90 GeV 
represents the region excluded after direct experimental searches [18]. 

As one can see, for tan (5 < 7 charged Higgs boxes give a positive correction to 
AMs, irrespective of the value assumed by the physical mass M^±. For this reason, 
we have chosen not to include such correction in the previous study on our mass 
scenarios: for each of them, Higgs contributions amount to a constant, positive shift 
of the final result. Then, in particular, the plots of Figs. 3-6 provide only a 'lower 
bound' on AM]^^, to be augmented by the contribution from Higgs boxes. 

It is interesting to give instead an indicative upper bound on the MFV MSSM 
contributions to AM^^ for low tan (3. To this aim, we can consider mass regime (A) 
- chargino dominated - which tends to give the largest positive corrections (around 
2 ps^^, see Fig. 3), together with Higgs contributions taken at a small value for 
Mjj± 100 GeV. By looking at Fig. 7, one can see that Higgs corrections amount to 

roughly another 2 ps~^. So, one can estimate MFV MSSM corrections to AM, not 
to exceed 4 ps~^. However, they are typically considerably smaller than this upper 
bound, as one can see by inspection of the various mass regimes (see also Fig. 8 
below) . This fact shows, for the case of Bs — Bs mixing, that explicit implementation 
of the MFV limit in the MSSM leads to a naturally small correction to the SM 
prediction: in fact, expansions (11) analytically realize the condition of "naturalness 
of near fiavour conservation" advocated in [12]. 

V. The case of AM^ 

We applied the entire strategy described above also to the case of AM^. Results are 
completely analogous, so we will limit to a few observations. 



-23- 




Figure 7: Charged Higgs boxes contribution to AMs as a function of the physical charged 
Higgs mass M^± . The different curves correspond to increasing values of tan /3 
from 3 (uppermost one) to 10 (lowermost one). The hatched area on the left of 
M^± = 90 GeV represents the experimentally excluded region [18] from direct 
search only. 



Keeping for the moment aside charged Higgs contributions, the magnitude of SUSY 
corrections AM^^ , normalized to the SM prediction AM^^, is basically the same 
as the corresponding quantity in the Bs case, in all our studied scenarios. Also the 
shapes of the distributions of values for the single contributions look very similar in 
the Bd and Bs cases. 

This leads to the following remark. Even in scenarios where gluinos give important 
contributions from scalar operators, the latter do not bring about a sensible depen- 
dence on the external quark flavours. In fact, leading scalar contributions go as ~ m^, 
and those in ~ mbUis or respectively ~ mfem^ are subleading ones. The effect of the 
latter is moreover completely hidden, in our case, by the lack of knowledge of the 
SUSY scales and of the MFV parameters. 

The inclusion of charged Higgs contributions does instead 'distinguish' the case of 
AMrf with respect to that of AMg- Now contributions to the LR scalar operator 
behave as ~ mi,mg(^d) tan^ (3 for AM^^^^) (look at the coefficient C^^^^ in the Ap- 
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pendix).^^ For tan/? < 10 , however, effects on tfie ratio AMj^/AM^ are witliin 1 %. 
Only for tan/3 > 10 do tfiey become larger than 3 %, and can be visible once the 
lattice error on ^ is at least halved with respect to the present value. 

6. Additional MonteCarlo's and role of tan (3 

In this Section we elaborate on our main findings, as described in the previous study. 
In particular, we provide further arguments in support of our results, by verifying 
them with a set of additional numerical studies. The latter are designed to explore 
possible loopholes in the above treatment, and to address the question how the picture 
changes when increasing tan (3 from strictly low values. 

From the above discussion, it is evident that the positiveness of the SUSY correction 
/S.Mf^ in the MFV MSSM is due to an interplay among different contributions, the 
most important being those from charginos and gluinos. In addition, for tan (3 < 7, 
Higgs boxes further shift the result by a positive amount, depending on the chosen 

Such findings followed from a MonteCarlo study in which SUSY masses were fixed 
and MFV parameters generated with flat distributions. To check for the robustness 
of our results, we also performed additional MonteCarlo's. 

• Random scan of both SUSY scales and MFV parameters 

In a first alternative set of runs, we scanned both SUSY scales and MFV parameters. 
In particular SUSY masses were assumed to follow flat distributions in the range 
Mi G [100,1000] GeV.^^ In such runs, we also set tan/3 to different values in the 
range tan/3 G [3,15]. This allowed us to study the resulting modifications in the 
various contributions, in particular in those from Higgs boxes, which for a light Mjj± 
tend to be negative when tan /3 > 7. 

The results of such global runs are reported in Fig. 8 for the cases of tan /3 = 3 
and 10. As the Figure shows, for tan /3 = 3 the positiveness of the sum of results 
is confirmed by the global run as well.^"^ The right panel of the first row also shows 
that, if M^± is assumed to be flatly distributed in the range M^± G [90, 1000] GeV, 
then Higgs contributions tend to give a distribution completely analogous to that for 
charginos (the latter hides almost entirely the former in the plot). 

The second row of Fig. 8 reports then the same distributions, when tan/3 = 10. 
As one can see, contributions from charginos, gluinos, neutralinos and mixed gluino- 
neutralino are left qualitatively unchanged by the variation in tan/3. On the other 
hand, Higgs contributions are very small and negative, as expected also from the 

^^At the same time, contributions to C^^^^ are suppressed by cot^ (3. 
^^The lower bounds for Mg and for M^± were set to 200 GeV and 90 GeV, respectively. 
^''We mention that a further verification was carried out by varying the 'B-parameters' for the 
effective operator matrix elements within the ranges allowed in [23]. Neither the positiveness of 

the sum of the results, nor the interplay bcwoon charginos and gluinos, resulting in the above 
regimes (A),...,(D), are touched at all by variations of the B-parameters. 
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Figure 8: Distributions of values for AMf^ in the MFV MSSM: sum of the contributions (left 
panels) and separate SUSY contributions (right panels). The distributions result 
from scanning the MFV parameters well as the SUSY mass scales (see 

text for details). First and second row plots refer to tan/3 = 3 and 10, respectively. 



lowermost curve in Fig. 7. However, this fact has quite a smah impact on the sum 
of all the results, which remains chiefly positive (see left panel). We mention that 
a further increase in tan/3 to 15 confirms a similar picture: all contributions but 
for Higgs boxes build up an almost totally positive distribution. The negativity of 
Higgs boxes, however, starts to matter, and the sum of all contributions is positive 
in roughly 86 % of the counts. 
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• Changing the ranges for MFV parameters 

A second set of runs was devised to check for possible variations in our findings, 
when the MFV parameters defining the expansions (11) are varied in ranges different 
from those specified in eq. (27). In the latter, the choice of 1 as the upper bound 
is dictated by various considerations. First, on 'aesthetical' grounds, such bound 
should be realistic, if the MFV expansion of the soft terms as functions of the Yukawa 
couplings is a 'natural' one. The interpretation is that and bi 'couplings' of 0(1) 
cause the new flavour violating effects originated by the soft terms to be at the 
same time CKM-like - i.e. generated by the SM Yukawa couplings, according to 
the very MFV paradigm - and natural [12, 1]. Second, on more technical grounds, 
enlarging too much the range of variation for the and bi may lead to a non-efficient 
exploration of the full parameter space allowed. However, on this latter point, we 
have verified that the shape of the distributions remains stable already after a few 
hundreds random values for the a^, 6,. 

In this second MonteCarlo study, we kept SUSY masses fixed to scenarios, as 
described in Section 5.1, and increased the upper bounds for the MFV parameters 
from 1 to 5. The main effect of increasing the ranges of variability is to 'smoothen' 
the difference between chargino and gluino distributions in the various scenarios and 
consequently to make regimes (A),..., (D) more similar to one another. However, the 
main features of each of them, as well as the positiveness of the sum of contributions, 
are left unchanged. 

On the whole, a choice of the MFV parameters according to the ranges (27) - 
mainly dictated by 'naturalness' - allows a better understanding of the interplay 
between chargino and gluino contributions in the different scenarios. It would be in- 
teresting to adopt a top-down approach to the determination of the MFV parameters 
ai,bi, by considering e.g. SUSY Grand Unified Theories (GUT) which at low energy 
reproduce the soft SUSY breaking structure of the MSSM. If such SUSY GUTs can 
be constructed to be minimal flavour violating [33], then the Oj, bi at the EW break- 
ing scale can be fitted by introducing the corresponding constraints (11) between soft 
terms and Yukawa couplings. The latter are both predicted within such models, in 
terms of the initial conditions at the GUT scale and of the running. 

We conclude this section by stressing that the feature of positiveness for the sum 
of SUSY contributions to AMg in the MFV MSSM is a precise signature of low 
tan P ^ 10. It is interesting to address the question whether a similar signature, but 
with sign reversed, applies to the MFV MSSM in the large tan/3 limit. This regime 
requires however consideration of an additional set of contributions, represented by 
the double neutral Higgs penguins [8, 34, 9, 10, 35]. This goes beyond the scope 
of the present paper. For an comprehensive related study in the context of grand 
unification, see Ref. [36]. 
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7. Various considerations on MFV 



The concrete application of the MFV hmit to a calculable NP model, like the MSSM, 
gives us now the opportunity to emphasize the main differences between the model 
independent approach of [1] and the former, phenomenological definition of MFV 
by [2]. We will also try and understand in which, among our studied scenarios for 
the SUSY scales, MFV contributions to meson oscillations are dominated by Qi, the 
operator already present within the SM. This case is referred to as 'constrained' MFV 
(CMFV) [5] and its study will allow us to understand how natural CMFV is within 
the MSSM. 

7.1. MFV: definition [1] versus [2] 

According to the definition of MFV by [1], all the flavour and CP violation at the 
EW scale is generated by the Yukawa couplings present in the SM. As we stressed 
several times above, this does not mean that new sources of flavour violation should 
be set to zero, but instead that they should be taken as functions of the SM Yukawa 
couplings. The functional dependence is in turn fixed by their formal transformation 
properties under the flavour group and gives rise to expansions like (11). 

On the other hand, the phenomenological definition of MFV by [2] does not start 
from the SM Yukawa couplings, but from the CKM matrix. In this approach, a theory 
is considered minimal flavour violating, when only interactions which display explicit 
proportionality to the CKM matrix are active. In this context, new sources of flavour 
violation are simply not considered, since they have a priori "nothing to do" with the 
CKM matrix. 

However, in the more general approach by [1], the CKM matrix becomes a by- 
product of the SM Yukawa couplings. Then the possibility of treating the Yukawa 
couplings as spurion fields, allows to classify in terms of them any new source of flavour 
violation, and this is why the definition by [1] has a completely general applicability. 
As a matter of fact, it also provides a framework for systematic routes to beyond- MFV 
(sec [37]). 

In the general context of [1], let us finally establish a contact with definition [2]. 
Since in the approach [2] new flavour violating structures are taken as unrelated to the 
SM Yukawa couplings, it is clear that one can reach this definition by setting the MFV 
parameters bi 0, in expansions (11). In this respect, an interesting MFV study 
where effects beyond the less general framework [2] are not visible has been reported 
in Ref. [38]. Specifically, in this case squark flavor mixing effects proportional to 
Yukawa couplings are negligible (in particular there are no gluino contributions) and 
the hi coefficients are consistently set to zero. 

7.2. Deviations from 'constrained' MFV 

In this Section, we finally study - in the context of the AB = 2 Hamiltonian - the spe- 
cial case in which MFV MSSM contributions are dominated by the operator Qi, the 
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one already present in the SM. This case corresponds to 'constrained' MFV (CMFV) 
[5]^^. We would like to stress that CMFV is a phenomenologically relevant limit: in 
fact, if MFV is motivated by the observation that experiments do not require new 
sources of flavour violation besides the SM Yukawa couplings, the constrained ver- 
sion of it is analogously justified by the fact that there is no compelling experimental 
evidence for contributions of effective operators other than the SM ones. 

The latter statement is actually strictly true when one supposes that new effective 
operators be "strongly coupled", i.e. multiplied by an overall 'coupling' factor taken 
to be 1 (or —1) [39, 1]. Bounds on new operator contributions are considerably 
alleviated, or removed altogether, when such operators arc weakly coupled, as is 
typically the case for SUSY contributions to the observablcs considered in the bounds. 

The limit of CMFV is a useful benchmark case to consider, both because of the 
generic argument on new operators mentioned above and because it is a very well sat- 
isfied feature of chargino contributions, which represent an important and sometimes 
dominating ingredient in the various mass regimes of MFV. In the present Section, 
we start again from our calculated MFV MSSM contributions to meson oscillations. 
We then compare - in the different mass scenarios considered - the part of the contri- 
butions which bears proportionality to Qi with the rest of the contributions, due to 
operators besides the SM one. This exercise will allow us to understand how 'natural' 
is CMFV within the MSSM. 

As a first step, we can define the ratio i?cMFV as 

_ lAMfP (Ci ^ 0)1 

= AiwF>d^ ' ^^^^ 

where Ci are the Wilson coefficients of the AB = 2 Hamiltonian (see eq. (23)) in the 
MFV MSSM. The ratio -Rcmfv provides a good 'quantitative' understanding of the 
relative importance of contributions from non-SM operators with respect to Qi, in 
the Bg case. 

Within every of the mass scenarios considered, the ratio i?cMFV is then a distribu- 
tion of values depending on the MFV parameters Oj, bi. In Figs. 9-12 we report such 
distribution in the four representative regimes studied in Section 5.2. In particular, 
left panels display i?cMFV for the sum of all contributions, leaving aside, again, those 
from charged Higgs, on which we will comment separately. The percentage reported 
on top of the plots gives the number of counts satisfying |i?cMFv| ^ 0.05. This bound 
can be considered - in our MFV MSSM case - a quantitative 'definition' of CMFV. 
According to it, a given estimate of H^^^ is taken to be CMFV, if contributions 
from operators other than Qi do not exceed 5 % of the contribution from Qi. We 
warn the reader that this upper bound (and the corresponding definition associated 
with it) has of course a large degree of arbitrariness: we chose 5 % in view of the still 
large hadronic uncertainties associated with the estimate of the {Qi}- 

As one can see from the left panels of Figs. 9-12, mass regimes (A) and (D) satisfy 



'For the sake of clarity, we mention that the definition of MFV adopted in [5] complies with [2]. 



-29- 




96.70 % CMFV 



-0.06 -0.04 -0.02 -0 ^0.02 0.04 0.06 0.08 0.1 




Figure 9: Ratio (29) of the contributions to AMf^ from operators other than Qi to that 
from Qi alone, for the set of SUSY scales chosen in Fig. 3, which is representative 
of mass regime (A). Left panel shows the sum of all contributions but for charged 
Higgs boxes, right panel reports the ratio for the separate contributions. 



well CMFV. Mass regimes (B) and in particular (C), on the other hand, do not.^^ 
Further insight can be gained by looking at the right panels of the same Figures, 
where Rqmfy is displayed for the single contributions (in this case the ratio (29) was 
calculated by restricting separately to Cf^ for gluino-gluino boxes, and so on for the 
other contributions). One immediately recognizes that regime (A) is Qi-dominated, 
since chargino contributions are the most important ones. Gluinos are not CMFV, but 
they are unimportant. In regime (D), instead, charginos and gluinos are of the same 
size, but both Qi-dominated. Finally, in regimes (B) and especially (C), CMFV does 
not apply: in fact, while charginos are always important and Qi-dominated, gluinos 
are not negligible and not Qi -dominated. As a matter of fact, Rcmfv for gluinos is 
in these cases very spread, typical values being in the range 20 -r- 40. We mention 
that such large values are not shown in the plots of Figs. 10 and 11. 

Let us then address how the above picture is modified when adding contributions 
from charged Higgses. The i?cMFV for the latter, in the case of tan /3 = 3, ranges 
between 0.035 and 0.015, for M^± = 90 GeV and 1000 GeV respectively. However, 
contributions from scalar operators steadily increase in importance when increasing 
tan/3: for tan/3 = 4 Higgs contributions are not CMFV if M^± ^ 500 GeV and, 
already from tan/3 = 5, in the full range of masses considered for M^±. Then, 
whether the total sum is CMFV or not, it depends on how large the contribution 
from Higgses is, with respect to the other contributions, i.e. on the Higgs mass 

^^We mention that the case displayed in Fig. 10 is the one with the highest CMFV percentage among 
the studied scenarios belonging to mass regime (B). 
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Figure 10: Same as Fig. 9, but for the set of SUSY scales, chosen in this case as in Fig. 4, 
which is representative of mass regime (B). 

chosen (see Fig. 7). 

We remark at this point that an increase in importance of scalar contributions 
affects, in general, AM^ and AM^ in a different way. However, in the case of Higgs 
boxes, we found that the ratio of AMd/ AMg has a sensitivity with respect to varia- 
tions of Mjj± below 1 %, when tan/3 is strictly small. Variations start to be visible 
only from tan/3 ^ 10, when they become > 3 %. This point was already stated at 
the end of Section 5, when discussing the case of AM^. 

Concerning contributions other than charged Higgses, we note that large deviations 
from Qi-dominated MFV apply when ^ is not small - regimes (B) and (C) However, 
such deviations turn out not to affect AMg and AM^ in a sensibly different way, so 
that their impact is not visible in the ratio AMd/AMg. The reason was, again, 
already explained at the end of Section 5. 

However, the above discussion allows us to place a warning message on the conven- 
tional wisdom that, for low tan/3, MFV in the AF = 2 Hamiltonian is Qi-dominated. 
This is clearly not true when Higgs and gluino contributions are included, and when 

is not small in magnitude. In this case, the approximation 1^ ^ 0, on the ground 
that tan /? is small, turns out not to work. 

7.3. A lower bound on AMg^a from CMFV 

Recently, it has been shown analytically in [40] that, within the CMFV models, 
the lower bounds on AM^^s are simply given by (AMd^s)sM- The proof holds, as- 
suming general exchange of charged gauge bosons, Goldstone bosons and physical 
scalars in boxes together with Dirac fermions, and under phenomenologically real- 
istic assumptions on their masses. On the other hand, two possible exceptions to 
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Figure 11: Same as Fig. 9, but for the set of SUSY scales, chosen in this case as in Fig. 5, 
which is representative of mass regime (C). 

the argument were found to arise in the presence of Majorana fermions and of U{1) 
neutral gauge bosons in box diagrams, which could individually bring AAI^^s below 
(AMrf^s)sM- Within the model-independent approach of [40], it seems impossible to 
exclude CMFV models involving negative contributions to AM^^s, caused by Ma- 
jorana fermions and/or U{1) neutral gauge bosons, that are not cancelled by the 
remaining positive contributions from the other possible particle exchanges. 

Within our study, we found in the previous Section that only mass regimes (A) and 
(D) are Qi-dominated, i.e. CMFV. In particular, the above mentioned exception, of 
a CMFV case in which Majorana fermion contributions are not negligible, arises in 
regime (D) (see Fig. 6). However in this case, Majorana contributions, represented 
by the gluino boxes, are also positive, so that, at least for low tan/3, the lower bounds 
of [40] still hold. Finally, regime (A) exactly corresponds to the assumptions made 
in the proof of [40], since Majorana contributions are negligible. In this case we 
numerically confirm the positivity of the NP contributions. 

We find actually remarkable that, within the MFV MSSM at low tan /3, the condi- 
tion AAIg^d > (AMs ,i)sM holds even when new operators give non-negligible contri- 
butions. 

8. MFV-Unitarity Triangle 

In this Section, we would like to emphasize that the so-called Universal Unitarity 
Triangle (UUT) [2] is valid only within the CMFV models, and not within the general 
class of MFV models, like the MSSM considered here. 

Indeed, the usual construction of the UUT is based on the value of the angle /?, 
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Figure 12: Same as Fig. 9, but for the set of SUSY scales, chosen in this case as in Fig. 6, 
which is representative of mass regime (D). 



measured by means of the S^Ks asymmetry, and on the value of the side Rt, obtained 
from the ratio AM^/ AMg, that within CMFV is independent of any new physics 
contributions. 

As already stressed in [10, 5], in the context of the MSSM at large tan/5, the 
presence of new operator contributions to AMg and AM^ generally modifies the 
relation between Rt and AM^/ AMg, so that it now reads [10, 5] 



Rt = 0.913 



\ ^ 1 


/l7.8/ps / 


1.23 


V AM, V 



AMd , l + /s 

^ R^^= , (30) 



0.507/ps ^ ' l + fd 



with AMg = (AMg)sM(l + fg) and ^ being a known non-perturbative parameter, 
C = 1.23(6). 

Basing on the results of the previous sections, formula (30) applies to the MSSM at 
low tan/3 as well, with Rgd 7^ 1, due to the presence of new operators that differently 
affect AM, and AMd- This can already be seen by adding to the SM expression for 
the mass differences, the contributions from the charged Higgs boxes, which do not 
depend on the MFV parameters ai,bi. As we showed in Section 7, when addressing 
deviations from CMFV, the ratio AMd/ AMg is subject to variations at the percent 
level, depending on the value chosen for M^± . Similar effects occur in principle when 
considering the other contributions as well, but in this case the variation is completely 
hidden by the lack of knowledge of the values for the MFV parameters Cj , bi and of 
course of the SUSY scales. The departure from CMFV is more pronounced in the 
MFV MSSM for large tan /5, but such analysis is beyond the scope of the present 
paper. 
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Figure 13: The blue area represents the correlation between Ri, and 7, valid in MFV. The 
orange area corresponds to the region with Rt fixed to the value in eq. (30). 
CMFV occurs in the intersection between the two areas, displayed in red. Finally, 
the green dashed box shows the la-allowed range for Rb and 7, as measured from 
tree-level decays. 

The above argument clearly demonstrates that the MFV MSSM, even at low tan (3, 
does not belong to the class of CMFV models, and consequently, the UUT obtained 
from Rt and f3^Ks is not generically valid in MFV. 

This discussion suggests that the unitarity triangle valid for all MFV models is 
not the UUT of [2] - and analyzed in detail by [4] - but a triangle constructed from 
the angle (5^Ks ^-^d the new-physics-independent value of \Vub\ or 7 from tree-level 
decays. 

In this context, the known value (i-t^Xs — (21.2 ib 1.0)° establishes a correlation 
between the side Rh, 

of the unitarity triangle and the angle 7, that is valid for all models with MFV. 
Such correlation is represented in Fig. 13 as a blue area, under the assumption 
/3 = l3^Ks = (21.2 lb 1.0)°. The orange area in the Figure displays instead the region 
characterized by Rt fixed to the value in eq. (30), with Rg^ = 1 and ^ = 1.23(6). 
The intersection between the two areas - displayed in red - is then the one allowed 
to CMFV. Fig. 13 also shows the Icr-allowed range from tree-level decays, namely 
62° < 7 < 102° and 0.40 < -Rfc < 0.45, as a green dashed box. The latter overlaps 
with the higher branch of the MFV area but not with the CMFV one. This indicates, 
on the one hand, somehow a 'tension' [5] between the tree-level determination of R^, 
and the one favoured by CMFV. On the other hand, the overlap between the green 
box and the blue area suggests that this tension disappears within MFV, provided 
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7 > 80°. 

It will be interesting to monitor with the help of Fig. 13 the progress in the 
determination of and of the angle 7 from tree-level decays, and to verify whether 
the correlation in question is satisfied by the data. 

Equivalently, with precise values of 7 and i?^, that are used for the construction 
of the reference unitarity triangle, it will be possible to find out whether the CMFV 
UUT or the MFV-UT is chosen by nature, or instead if one has to introduce non- 
MFV interactions to fit the data. This would occur if the experimental point (7, Rb) 
lies outside the blue area in Fig. 13. 

9. Conclusions and Outlook 

In this paper, we have applied the effective field theory definition of Minimal Flavour 
Violation (MFV) to the MSSM and explicitly shown how, by this definition, the new 
sources of flavour and CP violation present in the MSSM become functions of the 
SM Yukawa couplings. We have subsequently applied the MFV limit to the MSSM 
AB = 2 Hamiltonian at low tan /3. 

Our findings can be summarized in the following main messages: 

1. In the MFV limit of the MSSM, the soft breaking terms become functions of the 
SM Yukawa couplings. This non-trivial functional dependence causes flavour 
violation due to soft terms to be not zero in MFV, but instead 'CKM-like'. 
The constraints imposed on the soft terms by the MFV assumption lead to a 
significant increase in the predictive power of the model. 

2. The super symmetric corrections to AMg^a low tan /? are found to be always 
positive with respect to the SM formula. This feature is due to an interesting 
interplay between the chargino and gluino box diagram contributions. 

3. Even at low tan (3, the MSSM does not in general belong to the class of models 
with CMFV, in contrast to the statements made in the literature. The presence 
of gluino box diagram contributions necessarily brings in new operators beyond 
the {V — A) ^ {V — A) one, whose importance depends on the mass regime 
chosen. 

4. The side Rt, used in the determination of the Universal Unitarity Triangle, is 
actually not a good constant in MFV. Within the MFV MSSM at low tan /?, 
we have found variations that can reach the percent level, depending again on 
the mass regime chosen, and on the value of tan /3. To resolve such variations, 
one however needs a lattice determination for ^ with an error at most half the 
present value. In the case of the MFV MSSM at large tan (3 the situation is 
in this respect 'easier', since larger deviations from the CMFV value of Rt are 
more likely. 

5. A unitarity triangle valid in all MFV models (MFV-UT) can be constructed 
using only \Vub\ or 7, from tree-level decays, and the angle /3, extracted from 
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S^Ks- particular, with the measured value of (3^Ksj MFV implies a testable 
correlation between \Vub\ and 7. With the present high value of \Vuij\, MFV 
favours 7 > 80°. The LHC program on the measurement of 7 is then of utmost 
importance to cleanly consolidate - or disprove - the MFV paradigm. 
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A. Wilson coefficients for the AB = 2 effective Hamiltonian 
in the MSSM 

A.l. Contributions 

Below we list the non-vanishing new physics contributions to the Wilson coefficients 
for Bg mixing, eq. (23). The case of mixing is obtained by the replacement 2^1 
(or 5 — > 4 where applicable) in the external quark indices. In the following expressions, 
it is always understood that internal indices are summed over their respective ranges, 
i.e. /, J, K = 1, 3 for quarks, i,j = 1, 6 for squarks, a, 6 = 1, 2 for charginos and 
a,b = 1, 4 for neutralinos. 

Charged Higgs contributions 

Ci^"^ = l||2^^3if-/3i^7*2i^}2^^^^ {2M^Do{ml,ml, Ml^,M^) cot'' /3 
- D2 {ml^ , ml^ , M^± , ± ) cot^ /3 - ID; {m^ , m^^ , M^± , M^) cot^ /?} , 

4 2 2 

r<H+H+ 92 T^* T^* f^s'f^h 
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a 



H+H+ 



4 



167r2 " ^2 J2 gj^^4^ 



2L>o(m2^,m2^,M2±,M2,)} , 



9t 



W 



•}• 



a 



4M^ 



+ 



^o(m^^ , , M^± , M^) (tan^ /? + cot' p) ) | 



The above coefficients C^^^^ and C^*^^ arise from the sum of — H^, — H^ 
and — boxes (with H2 the charged Goldstone boson, entering calculations 
away from the unitary gauge). 



Chargino contributions 

Cf""" = -4^2D2{Ml,Mi^,Ml,Ml^ [V^^,]^^ [V^^,]^^^ [V^^^, 
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Cfx^ = --^^D,iM^^,Ml,^,Ml,Ml) [y,i,,]^^3 [V^,,l^^ [Vi^.r^^ 

^^'""^ = ^^2«'^^.'<'0 [VxU^]a^3 i^xU^^S [^xU^t,2 KuSm2 ' 

CX-x- ^ -M^DoiMl,,Ml,^,Ml,Ml) [vSraU Kuah.^ Kua^ ^IL^ 

The coefficients Cf ^ and ^ can be obtained from ^ and C3 respec- 
tively by the replacement of left- and right-handed chargino couplings V^. 
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Neutralino contributions 



= -^^2(m|,,,mB^,mJ,,mJ.) [vir>,]l,^ [f,%],^.3 [k,V]^^3 



- -^^^'^^"^^^^"^pKl^Kl) \^xD<i]u2 [^ioa],,, [Vina]^,, [V^Da] 



ai3 



= ^^^i^o(M|.^,M|,^.,M^o,M^o) , 



^5^°'° = -^^2(M|,^,M|,^,M^o,Mj) [y^^^j;^ [y/^j;.^ [y,^^,]^^3 [i^«,J^^.3 



- ^;^^o(M|.^,M|,^.,Mjo,Mjo) [y,^^,];., [V^U,,s KU 



ai3 



Interchanging left- with right-handed neutrahno couphngs in C12 3 yields the expres- 
sions for the Wilson coefficients C12 3 



Mixed neutralino and giuino contributions 



167r2 6 --^-i^v-iJj^-x-a 

[{ZD)3i{ZD)z3 [V^DdVa^ fed] + (^i?)2i (^d)2j fed] feJa^^} ' 



^|X0 ^ ^^Do(M|,,M|.,,M^o,M,) {3(Z.)3,(Z,)^, [l^,-,J„^3 

ai2 [^xDd]aj2 (-^0)51 (-^0)5 j [^;^d] aj3 [^xi'Jaj-3} > 
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- {ZD)3j{ZD)&i [K^bJla [^xDd]*aj2 ~ (■^^)2j (^b)!^ „i3 [^xDd\aj3 

- KZDhj{ZD)Qi [^xoJls [^X™]Ii2 ~ ^iZD)*2j{ZD)li [V^Dd] ^^3 ["^x^J 



To obtain the coefficients Cf^^ from Cf 2% one again has to interchange left- with 
right-handed neutrahno couphngs. In addition, analogous replacements have to be 
performed for the elements of the down squark mixing matrix that explicitly appear 
in the above expressions: {ZD)2i <-> {ZD)5i and {Zd)^^ (2'ip)6i. 

Gluino contributions 



= -y||2^^I^0(M|,^,M|,^.,M|,M|)(Zo)3^(^Z.)3,(^d)2.(^Z^)2, 

- ^^I?2(M|,^,M2^.,M|,M|)(Zi5)3i(^i?)3,(^D);,(^D)^,- , 



Cl' = -^^MlDo{Ml^,Ml,^,MlMl){ZDMZnhjiZD)UZD)^ 
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C~f = ^lMlDo{Ml^,Ml^,Ml,Ml){ZDMZDhjiZD)UZD)lj , 



^ 93 
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+Y||2Y^2(M|,^,M|,.,M|,M|)(ZB)3i(Zi,)6,{3(ZB);,(Zz,)^,-2(Z^ • 

To get the corresponding expressions for the coefficients 2 3 one again has to carry 
out the following replacements for the down squark mixing matrix: (.^0)21 ^ {^0)5% 
and {ZD)3i ^ {ZD)6i- 

A. 2. Chargino and Neutralino couplings 

Here we explicitly report - in terms of gauge couplings, Yukawa couplings and rotation 
matrices - the chargino and neutralino couplings used in the expressions for the 
Wilson coefficients. 

\yxUd\aiI = (^- — iZuTKiiZ+)la + {Zu)lK+3)iiZ+ho^ ^KI , 



[^xDdlail = ~3^*^'^''^)(^+3)»*^^^)la + [Z d) Ii{Z Ny3a 



A. 3. Loop functions 



Finally, we give the explicit expressions for the loop functions that appear in the 
Wilson coefficients listed above. 

n r 2 2 2 2n mjlnml 



{m\ — mf)(m3 — mf)(m2 — mf) 

+ {1 ^ 2} + {1 ^ 3} + {1 ^ 4} , 



D,{mlmlmlml) = \ 



mf In mf 



(777,4 ~ Tni){m^ — mi){m2 — mi) 

+ {1 ^ 2} + {1 ^ 3} + {1 ^ 4} 
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